Measurement of Muon Capture on the Proton to 1% Precision and 
Determination of the Pseudoscalar Coupling g p 
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The MuCap experiment at the Paul Scherrer Institute has measured the rate As of muon capture 
from the singlet state of the muonic hydrogen atom to a precision of 1 %. A muon beam was 
stopped in a time projection chamber filled with 10-bar, ultra-pure hydrogen gas. Cylindrical 
wire chambers and a segmented scintillator barrel detected electrons from muon decay. As is 
determined from the difference between the fi~ disappearance rate in hydrogen and the free muon 
decay rate. The result is based on the analysis of 1.2 x 10 10 fx~ decays, from which we extract 
the capture rate A s = (714.9 ± 5.4 sta t ± 5.1 sys t) s _1 and derive the proton's pseudoscalar coupling 
9p(lo = -0.88ml) = 8.06 ± 0.55. 

PACS numbers: 23.40.-s, 24.80. +y, 13.60.-r, 14.20.Dh, 11.40.Ha, 29.40.Gx 



We report a measurement of the rate Ag of ordinary 
muon capture (OMC), 

H~ + p n + Up, (1) 

from the singlet state of the muonic hydrogen atom. The 
analysis uses the complete dataset of the MuCap experi- 
ment, with significantly smaller systematic and statistical 
uncertainties compared to our earlier publication pQ. 

For the low momentum transfer q$ = —0.88m 2 in pro- 
cess ([I]), the Standard Model electroweak interaction re- 
duces to an effective Fermi interaction between the lep- 
tonic and hadronic weak currents. While the leptonic 
current has a simple 7^(1 — 75) structure, the hadronic 
current between nucleon states is modified by QCD, as 
expressed in a model-independent way by the introduc- 
tion of form factors. Since second-class currents are sup- 
pressed, muon capture on the proton involves <?y(<?o) and 
9m{%)i t ne vec tor and magnetic form factors in the vec- 
tor current, as well as g A (q$) and g P (qo), the axial and 
pseudoscalar form factors in the axial current [2H1]. The 
first three are well known and contribute only ~0.4% 
uncertainty to the determination of A5 [5]. Our mea- 
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surement of As determines g P = g P {q$), the least well 
known of these form factors. 

The pseudoscalar term in the axial nucleon current has 
played a significant role in the understanding of weak 
and strong interactions. Initial estimates were based on 
the concept of a partially conserved axial current, fol- 
lowed by the recognition of its deeper significance as a 
consequence of chiral symmetry and its spontaneous and 
explicit breaking [6j. These ideas were foundations for 
explaining the generation of hadronic masses and the de- 
velopment of chiral perturbation theory (ChPT), the ef- 
fective field theory of low-energy QCD. Based on well- 
known low-energy constants, 

^Theory = g ^ ± Q ^ ^ 

was derived within ChPT [H[7], with good convergence to 
two-loop order [8J. Though lattice QCD has advanced to 
unquenched calculations of g A and g P [fH [10] , the preci- 
sion of the ChPT prediction in Eqn.[2]remains unmatched 
and stands to be tested experimentally. 

Muon capture on hydrogen is the most direct means to 
determine g p . Such experiments are complicated by the 
fact that negative muons stopped in hydrogen form not 
only \ip atoms, but subsequently ppfi molecules where 
the capture rate differs significantly. Prior to MuCap, 
the most precise capture rate was measured in liquid hy- 
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FIG. 1: (color online) Upper panel: (a) Schematic sideview 
of a typical muon stop and decay electron. Lower panels: 
zoomed- in special event topologies (b) Rare large-angle fi — p 
scatter event with ~ 10~ 2 //i probability, c) Very rare de- 
layed capture on impurity in fiO — > N*v with ~ 3.4 x 10~ 6 //^ 
probability. The interpretation of the event displays will be 
described in the text. Images not to scale. 

drogen (LH2) where ppfi forms rapidly. The value 
of g P extracted under these conditions depends critically 
on the poorly known ortho-to-para ppfx transition rate, 
X op [121 [13] . The related, rare radiative muon capture 
(RMC) process fi~p — > nuj is less sensitive to A op , but 
the first experimental result on g p |14j disagreed with 
theory. For a discussion of this puzzling situation and 
muonic processes in hydrogen, see PHI]. 

The MuCap experiment was designed to significantly 
reduce the density-dependent formation of ppfi molecules 
by employing a gas density of <j> « 0.01 (relative to LH2). 
In these conditions ^97% of the muon captures occur 
in the fip singlet state. The experimental concept is 
sketched in Figure [lji. A 34 MeV/c muon beam was 
stopped in a time projection chamber (TPC) filled with 
10-bar, ultra-pure hydrogen gas [T5]. The TPC was used 
to discriminate between muons that stop in the gas and 
those that reach wall materials, where capture proceeds 
much faster than in hydrogen. Arriving muons were de- 
tected by an entrance scintillator (/iSC) and a propor- 
tional chamber (/xPC), and tracked in the TPC. Outgo- 
ing decay electrons were detected by concentric multi- 
wire proportional chambers (ePCl and ePC2) and a seg- 
mented scintillator barrel (eSC). The decay times were 
histogrammed and fit to an exponential. The difference 
between the observed disappearance rate A M - and the 
free muon decay rate A M + [16] is attributed to muon cap- 



ture, A s « A M - - A M + . 

The new data reported here were collected during the 
2006 (R06) and 2007 (R07) running periods in the ttE3 
muon channel at the Paul Scherrer Institute. Properties 
of the new datasets are compared to our published result 
(R04) PQ in Tabled Besides the 1.2 x 10 10 muon-electron 
pairs from /i~ stops in hydrogen, additional systematic 
data included 0.6 x 10 10 fi + decays and fi~ data collected 
when the target gas was doped with elemental impurities 
(nitrogen, water and argon). 

TABLE I: Main features of MuCap production runs. Statis- 
tics of fully reconstructed fi— e pairs, deuterium concentration 
c D , water concentration c H2 o determined by humidity sensor 
(not present in R04), and observed impurity capture yield per 
muon, Yz. 



quantity 


R04 


R06 


R07 


Statistics 


1.6 x 10 9 


5.5 x 10 9 


5.0 x 10 9 


c D (ppb) 


1440 


<60 


<10 


c h 2 o (PPb) 


n/a 


18. 


8.7 


Yz (ppm) 


12 


6.3 


3.4 



Although the experimental methodology closely fol- 
lowed that of our first result, several hardware upgrades 
implemented between R04 and R06 led to significantly 
enhanced performance. In R06 the TPC was operated 
with ~2.5x higher gas gain than in R07. As this affects 
critical chamber parameters, the comparison of the two 
runs provides an invaluable consistency check. 

Events with multiple muons in the TPC (pileup) need 
to be rejected as they distort the extracted disappearance 
rate. The maximum rate of the DC muon beam employed 
in R04 was throttled to minimize pileup. In R06 and 
R07, the loss of events to pileup was largely eliminated 
by the introduction of a 25-kV, fast- switching electro- 
static kicker [17] . The detection of a muon traversing the 
fiSC triggered the kicker, which deflected the beam for a 
period of 25.6 fis. The beam extinction factor was « 100 
and the rate of pilcup-free data was three times larger in 
R06/7 than in R04. 

Another essential change was the greatly improved iso- 
topic and chemical purity of the TPC gas. When deu- 
terium is present (concentration c D ), muons can form 
fid atoms, which, due to a Ramsauer-Townscnd mini- 
mum in the scattering cross section [3] , can diffuse out of 
the fiducial volume and distort the disappearance rate. 
To separate the hydrogen into its isotopic components, 
a new cryogenic distillation column was installed. Peri- 
odic gas samples were analyzed externally using atomic 
mass spectrometry [TS] . For the limits on c D listed in Ta- 
ble HJ transfer to fid leads to distortions of < 0.74 s _1 and 
< 0.12 s _1 for the R06 and R07 run periods, respectively. 
A higher sensitivity of the atomic mass spectrometer was 
responsible for the improved limit in R07. 

In the presence of Z > 1 impurities, muons prefer- 



3 



entially transfer from \ip to \xZ atoms, distorting the 
disappearance rate. Extended baking of the TPC and 
increased flux and filtering in the continuous gas purifi- 
cation system |19j led to a fourfold reduced impurity level 
compared to R04. Moreover, the installation of a humid- 
ity sensor before R06 allowed monitoring of the dominant 
chemical impurity. 

Our novel hydrogen TPC was key to the experiment. 
In the sensitive volume of i x y x z = 15 x 12 x 28 cm 3 
with an applied field of 2 kV/cm, ionization electrons 
from stopping muons drifted down towards the readout 
plane with a velocity v y — —5.5 mm/^s. They were am- 
plified in a multi-wire proportional chamber region with 
72 anodes perpendicular and 36 strip wires parallel to 
the beam. Anode and cathode signals were discriminated 
with three energy thresholds and read out by TDCs in 
200 ns time intervals. The event display in the center 
of Figure [IJi shows the y-z projection of a typical muon 
stop in the TPC. The threshold E t m 15 keV is indi- 
cated by green pixels. Blue pixels denote the threshold 
E 2 ~ 55 keV, which was set just below the muon's Bragg 
peak. The threshold E3 sa 315 keV (red pixels in Figure 
[lp) was set to record nuclear recoils from muon captures 
on impurities. In addition to the primary TDC-based 
readout of the TPC, new 12-bit, 25 MHz flash analog-to- 
digital converters recorded selectively triggered events. 

About 30 TB of raw data was processed at the 
NCSA supercomputing facility in a multi-stage proce- 
dure. Muon stop candidates were constructed from the 
/xSC time and TPC track information. Muon pileup 
events, flagged by the entrance counters, were rejected. 
Because the combined inefficiency of the entrance coun- 
ters was less than 10 -4 , the residual pileup distorted 
the observed muon disappearance rate A„- by less than 
0.5 s _1 . Contiguous pixel regions in the TPC were then 
fit to a straight line, as indicated by the black line along 
the muon trajectory in Fig. [T^,. In the case of large an- 
gle scattering (Fig. [TJd) a two-line fit was applied. The 
muon stop location (red circle) was identified as the most 
downstream E% pixel. The muon track requirements were 
optimized so as to minimize possible distortions to A M - 
while suppressing events where the muon could have left 
the hydrogen gas. Muons that stopped within a fiducial 
volume Ax x Ay x Az = 10.4 x 8.0 x 20.4 cm 3 were ac- 
cepted. The minimum track length was 3.2 cm, and the 
maximum fit \ 2 was 2. A„- was stable against variation 
of the track length and x 2 cuts. However, variations in 
the fiducial volume boundaries produced statistically dis- 
allowed deviations, for which a systematic uncertainty of 
3.0 s^ 1 was assigned. 

Electron tracks were constructed from coincidences be- 
tween an eSC segment (comprising four photomultiplier 
tubes) and hits in the two ePCs (each requiring an an- 
ode and at least one cathode plane). In the R06 and 
R07 run periods, the time and gain stabilities of the eSC 
were verified by recording their signals in 8-bit, 450 MHz 



waveform digitizers. While the TPC gain was insufficient 
to produce electron tracks with contiguous pixels, a vir- 
tual track in the TPC was reconstructed from hits in the 
eSC and ePCs, as indicated by the red line in Figure [ij,. 
A cut was placed on the impact parameter b between 
muon stop and electron vector. Tighter cuts on b produce 
lower backgrounds, but they introduce a time-dependent 
acceptance due to [ip diffusion. Although atoms dif- 
fuse only at the mm scale, changes in A„- vs. b were 
observed. This behavior fixed the single parameter of a 
/ip diffusion model in good agreement with theory [20) . 
For the loose applied cut of b < 120 mm, the model was 
used to determine small corrections (—3.1 ± 0.1) s _1 and 
(-3.0 ± 0.1) s" 1 for R06 and R07, respectively. To check 
that A^- was insensitive to the electron track definition, 
we also constructed coincidences requiring different com- 
binations of anode and cathode planes within the ePCs. 
This revealed slightly non-statistical variations in A„- , 
which were fully covered by a 1.8 s _1 systematic uncer- 
tainty. 

In extreme cases of example Fig. [TJd, muons scatter 
through large angles, leave the TPC volume, and stop on 
surrounding materials. Because of the lower TPC gain 
during R07, there were often gaps in the tracks of scat- 
tered muons, making it difficult to reliably identify these 
events. Moreover, the recoil proton could deposit enough 
energy in the TPC to trigger the E2 threshold, mimick- 
ing an acceptable muon stop. However, these events were 
unlikely to deposit enough energy at the scattering ver- 
tex to exceed the ^-threshold on neighboring anodes. In 
the analysis of the R06 and R07 datasets, we required at 
least two consecutive E 2 anodes at the end of the muon 
track. 

This cut introduced a subtle systematic effect. Elec- 
trons that traversed the muon's drifting ionization charge 
occasionally deposited enough additional energy to ele- 
vate a muon's E\ signal above the E2 threshold. In rare 
instances a muon stop with a single E2 anode was pro- 
moted to a stop with two neighboring E 2 anodes. Such 
events would pass the /1 — p scatter cut described above 
with a decay-time-dependent acceptance and therefore 
distort the extracted disappearance rate. Because posi- 
tive muons are sensitive to the charge interference effect 
but do not capture on nuclei, we were able to measure the 
induced distortion [31]. The method was supplemented 
by neutron data collected in 8 large liquid scintillator 
detectors: muons scattered into Z > 1 materials were 
found to yield copious neutrons from nuclear capture. 
The resulting corrections were (— 12.4 ± 3.22) s -1 and 
(-7.2 ± 1.25) s" 1 for R06 and R07, respectively. The 
correction was sensitive to the E 2 threshold, which we 
decreased in R07, suppressing the interference effect. 

As illustrated in Fig. [lp, nuclear capture on impuri- 
ties was identified by the presence of an E 3 threshold 
signal in the TPC. This allowed continuous in-situ mon- 
itoring of the yield Yz of these events. The average val- 
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ues for Yz over the three datasets are given in Table [T] 
and track well with the humidity sensor readings. To 
calibrate the necessary correction, special runs were con- 
ducted in which the hydrogen gas was doped with known 
amounts of nitrogen or water vapor. The changes in the 
disappearance rate and Yz were measured relative to the 
pure, undoped hydrogen data. Scaling by the observed 
Yz then determined the corrections for residual impuri- 
ties: (-7.80 ± 1.87) s- 1 and (-4.54 ± 0.93) s" 1 in R06 
and R07, respectively. 

To obtain the final muon disappearance rate, the muon 
stops and electron tracks were first sorted into muon- 
electron pairs. The decay time, t = t c §c ~ ^sc, wa s 
histogrammed and fit with the function N(t) = Nq ■ w ■ 
X p - • e f- c + B over the range 160 ns < t < 19000 ns. 
The bin width w was fixed at 80 ns, while Nq, B and A^- 
were free parameters. To avoid analysis bias, the exact 
clock frequency was hidden from the analyzers. After it 
was revealed, we obtained 

\- (R06) = 455857.3 ± 7.7 stat ± 5.1 syst s" 1 , (3) 
\ t - (R07) = 455853.1 ± 8.3 stat ± 3.9 syst s" 1 . (4) 

Because the fit x 2 /DOF = 1.2 ± 0.1 was slightly larger 
than expected for both R06 and R07, these values reflect 
inflated statistical uncertainties following the ^-factor 
prescription [5]. The 3-parameter fitting procedure was 
complemented by applying a full kinetics fit which in- 
cluded all atomic- and molecular-state effects as well as 
water and nitrogen impurities; the result was consistent 
within 0.2 s _1 . 

In order to check the consistency of our result, we ex- 
amined changes in A„- with respect to variations in data 
selection. The fit start and stop times were varied over 
a range of several microseconds and the parameters re- 
mained stable. Only statistical variations were observed 
when the data was sorted chronologically by run number. 
Since many of the subtle couplings between the muon and 
electron definitions are geometrical, the observed stabil- 
ity of the result with respect to azimuth was a critical 
cross-check. 

Table HTl summarizes the aforementioned corrections to 
our A„- result as well as the systematic uncertainties. 
Two additional corrections are required to correctly ex- 
press A„- as: 

A p - = (A M+ + AA W ) + A s + AA pPM . (5) 

Here AA W is a calculable Lip bound-state effect [22l [23] , 
while AA ppfJi accounts for the ~3% of muons that cap- 
ture from molecular states. The latter depends on X op 
and A PPM and is derived from fits to simulated data 
generated with the precise experimental conditions (gas 
density <j> — 0.0115 ± 0.0001, background level, impu- 
rity concentrations, and fit ranges). We used X op — 
(6.6 ±3.4) x 10 4 [1] and a newly determined value of 
X ppl ^ = (1.94 ± 0.06) x 10 6 s -1 [23], which was measured 



TABLE II: Applied corrections and systematic errors. 



Effect 


Corrections and uncertainties [s 1 ] 




R06 


R07 


Z > 1 impurities 


-7.8 ±1.87 


-4.54 ±0.93 


Li—p scatter removal 


— 12.4 ± 3.22 


—7.2 ± 1.25 


fip diffusion 


o i i n 1 n 


— o.U ± U.1U 


/id aimision 


i n ha 
± U. /4 


1 n 10 


Fiducial volume cut 


±3.00 


±3.00 


Entrance counter ineff. 


±0.50 


±0.50 


Electron track def. 


± 1.80 


± 1.80 


Total A„- corr. 


-23.30 ±5.20 


-14.74 ±3.88 


fip bound state: AA MP 


-12.3 ±0.00 


-12.3 ±0.00 


ppp states: AA PPM 


-17.73 ±1.87 


-17.72 ± 1.87 



by admixing 19.6 ± 1.1 ppm of argon to the TPC's hy- 
drogen. It agrees with the previous world average [4] 
but is three times more precise. Evaluating Eqn. [5] 
and including the updated positive muon decay rate of 
455170.05 ± 0.46 s _1 [Tj5], we determine the singlet 
capture rates: 

A S (R06) = 717.3 ± 7.73 stat ± 5.55 syst s" 1 , (6) 
A S (R07) = 713.1 ± 8.33 stat ± 4.34 syst s~\ (7) 

We also update slightly our previous publication [T] us- 
ing the latest values for A^+j A op , and X ppll , to obtain 
A S (R04) = 713.5 ± 12.5 stat ± 8.6 syst s" 1 . Accounting for 
correlated systematics among these three datasets, we 
report a final, combined result 

A MuCap = n4 g ± 5 4tat ± g s -!_ (g) 

This new result is in excellent agreement with recent 
theory [2"5H2T| . With the latest PDG values for V ud = 
0.97425 and ^(0) = 1.2701 5J, we derive from Ref. [17] 
an updated A| h ' = 712.7 ± 3.0 ± 3.0 s _1 , using their 
g P = 8.2. Here, the two uncertainties stem from lim- 
ited knowledge of g A and radiative corrections. Based 
on Ref. [37J we then determine g P from MuCap via 
A MuCa P = A Th' (1 + . 6 265A 5A -0.0108(^ uCap -^ h )) 2 . 
The equation allows for a recalculation of our result if 
5^(0) changes by Ag A in the future. We derive 

ff MuCap^2 = _Q_ 8 g m 2 ) = g g ± q 4 g ± q 2g) ^ 

where the two uncertainties arise from the error propa- 
gation of A^ uCa P anc j j\Th ^ respectively. 

Figure [2] illustrates the excellent agreement with the 
theoretical prediction, Eqn. [2] and highlights MuCap's 
reduced sensitivity to the molecular parameter X op . This 
answers the longstanding challenge of an unambiguous 
measurement of g P , generated by the mutual inconsis- 
tency of earlier experiments (OMC, RMC) and their 
strong sensitivity to X op . Corroborating values for g P 
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FIG. 2: Extracted values for g P as a function of the poorly 
known molecular transition rate X op [121 1131 I28| . In con- 
trast to earlier experiments (OMC [TT], RMC [H]), MuCap 
is rather insensitive to this parameter. 



are obtained in recent analyses [29J [30] of an earlier 0.3% 
measurement of muon capture on 3 He |31j . with uncer- 
tainties limited by theory. MuCap provides the most 
precise determination of g P in the theoretically clean /ip 
atom and verifies a fundamental prediction of low-energy 
QCD. 
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